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Abstract. Gravitational wave (GW) detectors operating on a long time range can be used for the study of space 
distribution of sources of GW bursts or to put strong upper limits on the GW signal of a wide class of source 
candidates. For this purpose we propose here a sidereal time analysis to analyze the output signal of GW detectors. 
Using the characteristics of some existing detectors, we demonstrate the capability of the sidereal time analysis 
to give a clear signature of different localizations of GW sources: the Galactic Center, the Galactic Plane, the 
Supergalactic plane, the Great Attractor. On the contrary, a homogeneous 3D-distribution of GW sources gives a 
signal without features. In this paper we consider tensor gravitational waves with randomly oriented polarization. 
We consider GW detectors at fixed positions on the Earth, and a fixed orientation of the antenna. 


1. Introduction 

In this paper, we base our study of detection of gravita¬ 
tional waves (GW) on a phenomenological approach as 
formulated by Finn (2001). We explore the reality that 
could emerge directly from observation, without consid¬ 
ering specific models of GW emission. As has happened 
frequently in the history of astronomy when a new ob¬ 
servational window opens, GW astronomy may lead to 
a revision of theoretical models and to the development 
of new ideas. Even upper limits on the energy of GW 
events are very informative because they put restrictions 
on the predicted characteristics of GW sources and then 
on the physics of GW themselves. Energetic GW bursts 
often arise in regions where the gravity effects are strong 
(supernova explosion, neutron star merging...). They thus 
carry information about fundamental gravity physics in 
an unexplored regime of strong field. 

Our phenomenological approach is also justified by the 
fact that standard calculations of GW signals from super¬ 
nova explosions and merging of relativistic compact mas¬ 
sive objects essentially depend on the details of the last 
relativistic stages of the gravitational collapse that are still 
poorly known (Thorne 1987,1997; Paczynsky 1999; Finn 
2001 ). 
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Several GW detectors of the third generation are op¬ 
erating with long duty times. Five resonant GW detectors 
of the International Gravitational Event Collaboration 
(IGEC) already cover 1332 days of observation (Astone 
2001), i.e. about 90% of the total time between 1997-2000. 
In 2002, the laser interferometric detectors TAMA300, 
GEO600 and LIGO started to operate, though with 
shorter duty cycles. 

Usual method of data analysis for GW detectors con¬ 
sists of searching for coincidences of events between two or 
more detectors (Prodi et al. 2000). Here we consider an¬ 
other method — sidereal time analysis (STA) —. that con¬ 
sists of recording all GW events as a function of sidereal 
time. This possibility was mentioned by Weber (1969) in 
his pioneering work. The Earth’s rotation, measured with 
respect to stars, gives a natural unit of time for ground- 
based GW detectors. The advantage of the STA method is 
that the noise events do not correlate with sidereal time, 
while true GW signals, coming from a certain distribu¬ 
tion of GW sources, is tightly related to definite sidereal 
time. This happens because, for a fixed distribution of GW 
sources, a given antenna pattern produces a fixed ampli¬ 
tude shape along the sidereal time. Hence, for a given 
detection threshold, GW events will be detected only for 
specific sidereal time intervals. 

Here, considering tensor waves first, we demonstrate 
the capability of the STA method to give a clear signa- 
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ture of some realistic distributions of GW sources. For this 
phenomenological approach we must choose a very general 
form of the GW burst amplitude and then a distribution 
of GW sources. In section 2 we explain how the amplitude 
of individual sources is chosen and how the considered 
distributions of sources are built. In section 3 we give, for 
each distribution, the relative rate of GW events along 
sidereal time, by considering a Weber-type bar-detector 
(EXPLORER) and an interferometric detector (VIRGO). 
In section 4 we analyze the results to show how to use our 
results for a practical comparison with real observations. 

Our main target is to prepare the tools for the concrete 
analysis of forthcoming GW signals. Our FORTRAN pro¬ 
gram is available on request. 


2. Description of normalized sources 

2.1. Amplitude of an individual GW burst emitter 

Let us recall briefly how the count of events is made. 
For more details see Baryshev and Paturel (hereafter, 
BP2001). 

Depending on the energy transformed into GW by a 
source, one will observe a change in the reference length l. 
This is measured by the observed amplitude h 0 b s = A l/l. 
This signal will be counted as a detected event only if h 0 b s 
is larger than the limiting amplitude hu m that defines the 
sensitivity of the detector. 

Many of the expected signals may be presented as 
the ’’unit pulse of GW burst” introduced by Arnaldi & 
Pizzella(1979). It is a GW burst in a sinusoidal wave with 
amplitude ho, frequency vq and duration r g . For tensor 
GW, the observed amplitude h 0 b s at a distance r, pro¬ 
duced by energy E gw , is (see e.g. Pizzella 1989): 


ft o6s = 1.4xHr 20 (±^)(^^)( ~ 9W 
r v 0 




Hence, a GW source at a fixed distance r, to a first ap¬ 
proximation, is characterized by three main observable pa¬ 
rameters E gw , vo and r g , for a wide class of radiation pro¬ 
cesses. 

In the present paper, we consider the normalized am¬ 
plitude that may be written as (r in Mpc): 


h obs \G(ot, d)| 
K ^ r 


( 2 ) 


The term G{a, S) is a geometrical factor that defines 
the position of the source, with respect to the detector. 
The expression of this term is given for bar and interfer¬ 
ometric detectors in BP2001. The quantity ft* defines the 
type of source. For galactic weak sources we will adopt 
ft* = 10~ 2 ° or ft* = 10~ 23 . For powerful extragalactic 
sources of GW, like supernovae, we will adopt ft* = 10” 17 
or ft* = 10” 19 . Obviously, our calculations may be rescaled 
by using any other combinations of main GW parameters. 

Because the polarization angle cannot be predicted, it 
has been taken at random. The net effect of this unknown 


polarization is that the number of events, at a certain 
level above the noise, is simply reduced (for more detail 
see BP2001). This does not change the shape of the count 
rate with sidereal time. Note that the polarization makes 
it that one cannot infer the true amplitude of a GW event 
from the observed amplitude. Nevertheless, a large ampli¬ 
tude means both a large true amplitude and a favorable 
polarization. 

2.2. Distribution of sources 

We now define a group of individual GW sources. The 
sources inside such a group have a well-defined distri¬ 
bution. For instance, a group can be the center of our 
Galaxy. In this case each individual source is a galactic 
emitter (like neutron star mergers or unknown sources of 
GW pulses) distributed in a sphere representing the re¬ 
gion around the center of our Galaxy. Another example 
of a group can be a supercluster of galaxies. In this case 
each individual source is a galaxy. The group can be also 
a flat distribution of emitters with a given thickness and 
a given size. An example is de Vaucouleurs’s Local Super 
Cluster (de Vaucouleurs, 1956, 1975). 

Each group of emitters is characterized by: 

— Number of individual emitters. Generally, this number 
is arbitrary but, when one considers several groups, 
their relative populations must be fixed. 

— The rate of GW burst emission, e.g. in solar masses per 
unit of time. Because we calculate the relative number 
of detected events (in percent) for an arbitrary number 
of sources, the rate of GW emission for each individual 
emitter may be set to one. This assumes that each in¬ 
dividual emitter has the same probability of producing 
a GW event. 

— The coordinates of each individual emitter. We 
adopted equatorial coordinates. So, one can use an ac¬ 
tual catalog of emitters (like a catalog of real galax¬ 
ies) or a simulated catalog, using coordinates obtained 
from an artificial distribution of emitters. 

— Distance of each emitter. For a simulated sample the 
distance is given by the considered distribution. For 
an observed catalog of galaxies, the distance can be 
deduced from the radial velocity using a given Hubble 
constant. 

There are several natural space distributions of GW 
sources that may be directly tested by the STA method. 
Here we consider the following models of GW source dis¬ 
tributions that will be tested in the next section: 

1. Homogeneous 3D-distribution. We consider a homo¬ 
geneous distribution of 10000 emitters in a sphere of 
radius R = 10 kpc centered on the detector. Each emit¬ 
ter is considered as a galactic object with a reference 
amplitude ft* = 10 _2 ° or ft* = 10 -22 for EXPLORER 
or VIRGO, respectively. This distribution will be used 
to show that, on average, there is a constant rate of 
GW events along sidereal time. Obviously, the result 
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is the same when a homogeneous 3D-distribution of ex- 
tragalactic objects is used (we made this verification). 

2. Galactic center. The galactic center is represented by 
one group of 10000 emitters located at the 1950 equa¬ 
torial coordinates: 5174224.0 — 285500. The number 
of sources is assumed to decrease following a Gaussian 
distribution with standard deviation (radius) of er = 
O.OOIMpc. The adopted distance between the Sun and 
the galactic center is O.OlOMpc. The kind of emitter 
is not explicitly identified, but one can imagine that 
the galactic center could be the source of many strong 
GW events. A reference amplitude h* = 10” 2 ° is used 
for the kind of emitters detectable with EXPLORER 
and h* = 10 ” 23 for emitters detectable with VIRGO. 

3. The galactic plane. The distribution will be centered 
on the galactic center 5174224.0 — 285500 at a distance 
of O.OlOMpc from the Sun. The radius of the distribu¬ 
tion (radius of our Galaxy) is taken as 0.0125-Mpc. 
The adopted thickness of the disk is 2 kpc. We con¬ 
sider 10000 emitters with a reference amplitude of 
h * = 10” 2 ° or K = 10 ” 23 for EXPLORER and 
VIRGO respectively. Each emitter will be distributed 
homogeneously in this thin disk. Note that, because 
our Sun is not at the center of this distribution, the dis¬ 
tribution along the galactic longitude is not uniform. 
Again, for this distribution, we do not use a priori in¬ 
formation on the types of GW sources. 

4. The supergalactic plane. The dominant structure 
around the Local Group is de Vaucouleurs Local Super 
Cluster (LSC). The Virgo galaxy cluster is the domi¬ 
nant structure of the LSC. We will consider a homoge¬ 
neous distribution of 10000 emitters in a disk of radius 
25Mpc and height 5 Mpc, centered on the Virgo clus¬ 
ter itself. The north pole of the supergalactic plane is 
at l = 47.37 degrees, b = +6.32 degrees. The refer¬ 
ence amplitude of each emitter will be h* = 10” 17 or 

= 10 ” 19 for EXPLORER and VIRGO respectively. 
According to recent work by Norris (2002) the faint, 
slowly-cooling Gamma Ray Bursts could be localized 
in the supergalactic plane. Although the effect is not 
yet confirmed it is interesting to predict what would be 
the rate of sources localized in the supergalactic plane. 

5. Virgo and the possible Great Attractor. This distribu¬ 
tion is made of 1000 emitters (galaxies) in the Virgo 
cluster, centered at the 1950-equatorial coordinates 
5122800.0+124000. We assumed a distance of 17 Mpc 
and a radius of 1.7 Mpc. We expect that Virgo will pro¬ 
duce the most easily detectable events (this is the rea¬ 
son why the Franco-Italian detector is called VIRGO). 
However, it has been suspected that a very huge struc¬ 
ture, the Great-Attractor (Dressier et ah, 1987; here¬ 
after, GA), could be present at three times the distance 
of Virgo, but hidden by our Galaxy. From the acceler¬ 
ation induced on Virgo, one can assume that the GA 
has 27 more galaxies than Virgo. We thus added to 
our synthetic sample 27000 galaxies, centered at the 
1950-equatorial coordinates 5150000.0 — 600000. We 
assumed a distance of 51 Mpc and a radius of 5.1 Mpc 


Z 



Fig. 1 . Definition of the coordinates. The random X, Y, Z 
coordinates define the position of an individual emitter 
G (e.g., galaxy). The X-axis gives the origin of Right 
Ascension a (7 direction). The XY-plane defines the equa¬ 
torial plane used to define the declination S. Each ran¬ 
dom emitter G is accepted (or rejected) for obeying the 
searched distribution. 

For each emitter of both Virgo and the GA we as¬ 
sumed that the reference amplitude is h* = 10” 17 or 
hi, = 10 ~ 19 for EXPLORER and VIRGO respectively. 
With a high sensitivity detector, we expect that the 
GA could become the most efficient source of GW 
events. This could give a clue for the detection of this 
putative dense region. 

2.3. Practical method for constructing the samples 

The Cartesian X, Y, Z coordinates of each emitter are cho¬ 
sen at random, for a uniform distribution, or following 
a Gaussian distribution around the center, for a discrete 
group like a cluster. These coordinates are chosen accord¬ 
ing to Figure [l]. The X-axis gives the direction of the ori¬ 
gin of Right Ascension (7 point). The North equatorial 
pole is defined by the /Taxis. The X, Y , Z coordinates are 
then transformed into equatorial coordinates and distance. 
Then the coordinates are transformed into galactic or su¬ 
pergalactic coordinates, depending on the distribution we 
are constructing. An individual emitter is then accepted, 
or not, depending on its position with respect to the dis¬ 
tribution. This operation is repeated until the number of 
emitters has reached the desired value. 

We considered the following detectors: 

Detector Latitude L Azimuth <E > 0 Type 

EXPLORER 46.5 deg -39 deg ba^ 

VIRGO_43.6 deg-19 deg interf. 
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The considered sensitivity range for EXPLORER is 
10 -18 to 10~ 20 and for VIRGO 10~ 20 to 10 -22 The main 
useful characteristics are: the latitude on the Earth (let us 
recall that the longitude on the Earth does not intervene 
because the count rate of GW events is registered as a 
function of sidereal time), the azimuth of the main axis 
(the main axis of the bar or the x-arm for an interfero¬ 
metric detector) and the type of detector. 

3. Results 

We calculated the GW normalized count rates as a func¬ 
tion of sidereal time for the two detectors listed above. 
For each of the studied distributions we give a figure with 
three panels: 

1. Panel 1 (top): A Flamsteed equal-area projection that 
shows the distribution in supergalactic coordinates. 

2. Panel 2 (middle): The shape of the normalized count 
rate predicted for the EXPLORER bar detector along 
the sidereal time axis. We consider five sensitivity lev¬ 
els: lO" 18 ' 0 , 1CT 18 - 5 , 10" 19 ' 0 , 10- 19 - 5 and lO" 20 0 . The 
corresponding curves appear from the bottom to the 
top; the higher the sensitivity the higher the rate of 
events. 

3. Panel 3 (bottom): The shape of the normalized count 
rate predicted for the VIRGO interferometric detec¬ 
tor along the sidereal time axis. We consider five sen¬ 
sitivity levels: lO" 20 - 0 , lCT 20 ' 5 , lO" 210 , lCT 21 ' 5 and 
;U)-22.o. The corresponding curves can be identified as 
said for panel 2. 

The figures || to || give the results for respectively: 

— A 3D homogeneous distribution of sources within 10 
kpc, 

— Sources localized in the Galactic center, 

— Homogeneous distribution in the galactic plane, 

— Homogeneous distribution in the supergalactic plane, 

— Sources localized in the Virgo cluster and in the Great- 
Attractor. 

Each curve gives the relative count number of GW 
events for equivalent sidereal time intervals centered at dif¬ 
ferent sidereal times. The Y-axis is then a relative number 
of GW events per unit of time, i.e. a relative rate of events 
(in percents). We give an example.- the middle panel of 
Fig. s It gives the rate of GW events for a homogeneous 
distribution of GW emitters in the galactic plane, as it 
would be seen from the EXPLORER bar detector. From 
the intermediate sensitivity (ICY 18 ' 5 ) one can see that the 
rate of events detected around the sidereal time ST = 5 h 
(80%) is about twice the rate expected around ST = 21 h 
(35%). The corresponding points are marked with an open 
circle in Fig. By comparing these kinds of curves with 
observed ones, it should be possible to have a signature of 
the localization of the GW sources. Then, from the sen¬ 
sitivity of the detector, one can place some limits on the 
radiated GW energy. 


Supergalactic coordinates 




sideral time (hours) 



sideral time (hours) 

Fig. 2. 3D homogeneous distribution of sources within 10 
kpc. See section 3 for explanations. 
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Supergalactic coordinates 
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Fig. 3. Sources localized in the Galactic center. See sec¬ 
tion 3 for explanations. 


Supergalactic coordinates 



Supergalactic Longitude (degrees) 



0 5 10 15 20 


sideral time (hours) 



0 5 10 15 20 


sideral time (hours) 

Fig. 4. Homogeneous distribution of sources in the galac¬ 
tic plane. See section 3 for explanations. 
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Supergalactic coordinates 


Supergalactic coordinates 
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Fig. 5. Homogeneous distribution of sources in the super- Fig. 6. Sources localized in the Virgo cluster and in the 
galactic plane. See section 3 for explanations. Great-Attractor. See section 3 for explanations. 
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4. Discussion and conclusions 


References 


We have presented a method to analyze the output signal 
produced by GW detectors operating over a long time. 
This ’Sidereal Time Analysis’ method consists of record¬ 
ing all GW events as a function of sidereal time, in order 
to calculate their rate (number of events per unit of time) 
as a function of phase, from Oh to 24h, for different sensi¬ 
tivity thresholds. If one admits that the noise is uniformly 
distributed along the sidereal time, the statistics of ob¬ 
served events can be performed even for very low signal 
to noise ratios (less than one). 

A random distribution gives a flat count rate (Fig. |j). 
So, any sources distributed in such a homogeneous manner 
and close enough to be detected will appear as additional 
noise. Fortunately, the closest sources would have a galac¬ 
tic origin and would not lead to a flat count rate as shown 
by Figs. |] and [|. One can note the similarities between 
these two figures. The main characteristic of sources of 
tensor waves localized in our Galaxy (Galaxy center or 
Galaxy disk) is that their count rate shows a deep depres¬ 
sion at 21 h sidereal time for the EXPLORER bar-detector 
and at 15h and 22h sidereal time for the VIRGO interfer¬ 
ometric detector. 

In figure we see that the shape of the count rate for 
sources localized in the supergalactic plane shows almost 
opposite phase in comparison with Figs. and [|. This 
is due to the fact that the supergalactic plane is almost 
perpendicular to the galactic plane. This means that one 
can characterize the origin of a signal simply from the 
shape of its count rate. 

Finally, the last calculations for the Virgo cluster and 
the putative GA also need comment (Fig|j|). We consider 
first the EXPLORER detector (middle panel). It shows 
that, at low sensitivity (lower curves), one will detect the 
Virgo cluster first, with a shape similar to the one found 
for the supergalactic plane (i.e., one peak at lOh sidereal 
time and one peak at 23h sidereal time) and that progres¬ 
sively one will detect a dominant region (between 2h and 
12h with a big depression at 18h) when the GA dominates. 
So, when the sensitivity of the detectors are improved, it 
will be possible to test the existence of the GA. 

For the VIRGO detector (Fig](] bottom panel) the ef¬ 
fect is similar. We see a depression at 5h when Virgo dom¬ 
inates, and a depression at 13h when the GA dominates. 

The next step consists of using real signals to start such 
an analysis. In particular, it will be possible to investigate 
the important case of coincident detections from multiple 
detectors and the case where one has only some events, 
the critical question being at what number of events can 
one have a significant result about the distribution of GW 
sources. 
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